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CONSIDERATIONS  FOR  THE  DESIGN  OF  GROUND 
CLUTTER  CANCELERS  OF  WEATHER  RADAR 


Dusan  S.  Zrnic',  Said  Hamidi,  and  Allen  Zahrai 
1 .  Introduction 

This  final  report  relies  heavily  on  the  preceding  interim  report  and  is,  in 
fact,  its  continuation  and  update.  There  is  very  little  repetition  and  we  often 
reference  sections  and  equations  from  the  interim  report.  Four  main  topics  not 
covered  in  the  interim  report  are  dealt  with  herein. 

The  interim  report  studied  the  ground  clutter  properties  at  the  Norman  sites, 

and  several  canceling  schemes  were  investigated.  It  was  shown  that  clutter  cross 

2  2 

sections  40  to  50  dB  below  1  m  /m  are  typical  and  that  clutter  spectrum  widths 
range  between  0.1  and  0.5  m*s~^  with  the  mean  value  of  0.25  m-s"^  at  an  antenna 
rotation  rate  of  10°s”\ 

Design  goals  for  the  filter  were  50  dB  rejection  in  the  stopband  with  a 
narrow  variable  notch  width,  a  total  annihilation  at  zero  Doppler,  and  one  dB 
ripple  in  the  passband.  Such  goals  can  be  achieved  easily  on  contiguous  pulse 
trains  with  either  a  recursive  filter  or  a  finite  impulse  response  filter.  How¬ 
ever,  if  an  interlaced  pulse  transmission  scheme  is  utilized,  canceling  becomes 
difficult.  The  interlaced  scheme  is  useful  because  it  enables  unambiguous  measure¬ 
ment  of  reflectivity  (at  one  pulse  repetition  frequency)  and  from  it  a  correct 
assignment  of  ranges  to  velocity  estimates  that  are  obtained  at  another  pulse  repe¬ 
tition  frequency.  It  was  demonstrated  that  recursive  filter  with  initialization 
does  meet  design  specifications  when  the  pulse  train  contains  about  64  samples 
corresponding  to  a  dwell  time  of  about  50  ms.  Under  these  conditions,  similar 
performance  can  be  obtained  by  properly  weighting  time  series  data  before  a  Fourier 
transform  is  applied. 

Section  2  treats  problems  that  clutter  causes  in  the  second  trip  area  (ring) 
for  velocity  estimation.  These  are  increased  area  obscured  by  clutter  and  a  dis¬ 
advantage  of  second  trip  echo  powers  with  respect  to  the  first  trip  clutter  power. 

Assignment  of  correct  ranges  to  velocity  estimates  when  different  clutter 
filters  are  used  in  the  reflectivity  and  velocity  channels  is  treated  in  Section  3. 

The  thrust  of  this  report  is  in  the  investigation  of  the  recursive  filter 
performance  when  a  small  number  of  pulses  (eight)  is  passed  through  the  filter. 

This  is  treated  in  Section  4.  As  in  the  interim  report,  the  quality  of  filters 
performance  is  judged  from  the  biases  in  the  mean  velocity  and  spectrum  width 


estimates  obtained  by  a  pulse  pair  processing  algorithm  behind  the  filter. 

Finally,  in  Section  5  are  presented  some  engineering  design  considerations 
for  the  recursive  filter.  This  filter  is  presently  being  built  in  support  of  the 
NEXRAD  requirements. 

2.  Effects  of  Ground  Clutter  Ring  in  the  Second  Trip  Area 
Ground  clutter  ring  in  the  second  trip  area  impedes  velocity  estimation 

because  clutter  there  competes  with  weather  signals  that  have  a  range  squared 

disadvantage.  Furthermore,  the  area  covered  by  clutter  is  larger  in  the  second 

trip  than  in  the  first.  We  illustrate  this  by  way  of  an  example  on  NSSL's  Norman 

Doppler  radar.  The  first  trip  area  at  three  elevation  angles  for  which  clutter  is 

larger  than  -10  dBZ  is  shown  on  the  histograms  (Figure  1)  and  in  Table  1.  Also 

listed  in  the  table  is  the  extent  of  clutter  in  range  r  .  The  area  A2  covered  in 

the  second  trip  (for  an  unambiguous  range  r  =  115  km)  is 

d 

A2  =  -c2  +  2irrcra  (1) 

which  translates  to  over  10  times  as  much  as  in  the  first  trip  (see  Table  1)! 

The  percentage  of  the  obscured  first  trip  area  is  only  between  2.4  and  4.7  percent, 
but  the  percentage  of  the  second  trip  area  (3  r  ir)  that  is  obscured  ranges 

a 

from  11  to  16  percent.  Figures  2  and  3  illustrate  this  effect  for  the  elevation 
angle  of  0.8°. 

Range  advantage  of  close  targets  for  an  unambiguous  range  of  r  =  115  km  is 

a 

drawn  on  Figure  4.  In  order  to  have  accurate  mean  velocity  estimates  (bias  less 
than  1  m*s_1)  from  the  pulse  pair  algorithm,  the  weather  signal  must  be  10  dB 
larger  than  any  spurious  components  including  clutter  (Zrnic1  and  Hamidi,  1981). 
Consequently,  a  patch  of  clutter  at  a  range  10  to  15  km  with  an  equivalent  re¬ 
flectivity  factor  of  10  dBZ  competes  with  weather  signals  up  to  40  dBZ  in  strength 
when  these  are  at  125  to  130  km. 

3.  Range  De-Aliasing  of  Velocity  Estimates  when  Clutter  Filtering  is  Employed 
We  shall  now  point  out  some  difficulties  that  are  inherent  to  a  dual 

sampling  mode  if  one  tries  both  to  eliminate  the  clutter  and  to  assign  correct 
ranges  to  velocity  estimates.  To  illustrate  the  problem,  consider  a  simplified 
interlaced  scheme  that  uses  a  pulse  repetition  time  T$  for  velocity  estimation  and 
2T$  for  reflectivity  (Figure  5). 

This  scheme  is  conceptually  identical  to  the  one  used  on  NSSL  radars  except 
for  the  2Tg  reflectivity  PRT  (NSSL  radars  use  4Tg).  However,  by  restricting  our¬ 
selves  to  first  and  second  trip  echoes  only,  we  avoid  needless  complications  while 
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Figure  5.  Interlaced  pulse  transmission  scheme. 
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retaining  essential  aspects  of  the  problem. 

Figure  6  illustrates  how  the  present  range-unfolding  algorithm  works  in 
presence  of  ground  clutter  and  first  and  second  trip  storms  (Hennington,  1981).  ; 

The  algorithm  examines  three  possibilities  by  comparing  powers  in  the  reflectivity  T 

channel  for  each  resolution  volume:  1)  the  first  trip  power  is  larger  than 
the  second  trip  power  P ^  by  10  dB.  Then  first  trip  velocities  are  unbiased,  and 
any  overlaid  second  trip  echo  is  not  usable.  2)  P ^  -  10  P^  the  second  trip 
velocities  are  valid,  and  weakly  scattering  first  trip  targets  are  obscured.  r 

3)  The  powers  are  within  10  dB  from  each  other,  and  velocities  cannot  be  assigned  i 

to  targets  in  either  trip.  Radar  equipped  with  phase  diversity  could  recover  part  j 

of  this  area  in  both  trips.  The  recoverable  area  depends  on  the  signal-to-noise 
degradation  that  one  is  able  to  tolerate.  For  instance,  if  the  first  and  second  t 

trip  signals  have  equal  powers  and  one  is  made  incoherent,  the  signal  to  "white 
noise"  ratio  over  the  Nyquist  interval  becomes  0  dB. 

Introduction  of  a  ground  clutter  canceler  in  a  Doppler  channel  changes  the 
situation  because  the  powers  now  from  the  first  trip  clutter  are  different  in  the  f 

two  channels.  In  principle,  the  optimum  scheme  is  a  dual  frequency  radar  with 

identical  clutter  filters  in  the  coherent  reflectivity  and  velocity  channels  so 
that  the  power  comparison  algorithm  would  work.  This  is  so  because  with  two 
frequencies  the  dwell  time  for  reflectivity  estimation  can  be  made  long  enough  to  ► 

insure  satisfactory  filtering.  Thus,  the  shapes  of  filter  notches  and  ripples  in 
the  passband  can  be  made  almost  equal.  But  then  a  second  coherent  receiver  is 
needed  and  possibly  a  log  channel  for  A6C  (if  instantaneous  AGC  is  not  contem¬ 
plated).  Although  the  notch  and  ripple  can  be  made  equal,  the  two  filters  •' 

necessarily  have  different  Nyquist  velocities,  with  the  one  for  reflectivity  being 
smaller  (by  a  factor  of  2  for  the  example  on  Figure  5).  Consequently,  there  may 
be  powers  from  storms  with  velocities  that  have  aliased  into  the  notch  of  the 
filter  in  the  reflectivity  channel.  These  would  be  canceled  while  no  similar  i 

canceling  would  occur  in  the  velocity  channel.  In  such  cases,  the  algorithm  based 
on  power  comparisons  in  the  reflectivity  channel  would  be  erroneous. 

In  the  remainder  of  this  section  we  consider  range  assignment  of  velocity 
estimates  when  a  single  frequency  radar  is  used  in  an  interlaced  mode.  We  assume  * 

that  several  batches  of  higher  PRT  pulses  in  the  velocity  channel  are  interlaced 
with  short  bursts,  one  or  at  most  two  pulses,  in  the  reflectivity  channel.  A 
scheme  with  eight  pulses  for  velocity  (see  section  4)  and  two  for  reflectivity  is 
workable  with  a  possibility  of  crude  (difference)  canceling  in  the  reflectivity  » 

channel.  In  the  sequel  we  imply  that  the  reflectivity  channel  is  used  to  set  AGC 
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Figure  6.  Illustration  of  obscuration  in  the.  velocity  and  reflectivity 
channels  for  the  transmission  scheme  as  on  Figure  5. 


and  de-alias  range  overlaid  velocities  by  comparing  powers  in  the  two  channels. 

The  reflectivity  measurements  should  also  be  made  from  the  velocity  channel  except 
at  times  when  storms  are  overlaid.  Even  in  the  absence  of  second  trip  echoes,  the 
powers  in  the  two  channels  are  not  going  to  be  exactly  equal  because  they  are  not 
estimated  from  the  same  samples  in  the  interlaced  mode  (this  difference  is  statis¬ 
tical).  Moreover,  the  quantization  noise  and  slight  nonlinearities  between  the 
transfer  functions  of  the  two  channels  cannot  be  perfectly  accounted  for  (this 
difference  is  systematic).  To  make  this  distinction,  the  clutter  and  signal  powers 
from  the  reflectivity  channel  are  primed  (Table  2).  The  input  powers  in  the 

A 

reflectivity  channel  are  clutter  and  first  trip  echo  Pz^=P^+  p^  or  power  P ^  which 
corresponds  to  the  second  trip  echo  for  the  velocity  channel  (Figure  7).  In  the 
velocity  channel,  powers  in  general  are  clutter  from  first  and  signals  from  both 
trips  p/Pg+fyiy 

At  the  output  of  the  clutter  filter  in  the  velocity  channel,  the  clutter 
residue  power  is  Pr  and  in  the  reflectivity  channel  clutter  residue  is  P^. 

Table  2  illustrates  the  powers  when  first  and  second  trip  targets  are  sampled 
in  the  reflectivity  channel.  Both  trip  targets  are  sampled  simultaneously  in  the 
velocity  channel. 

Now  it  suffices  to  compare  the  power  in  the  reflectivity  channel  with  the 
filtered  power  Pyf  in  the  velocity  channel.  We  assume  that  all  the  powers  are 
sufficiently  higher  than  the  receiver  noise.  Thus  range  de-aliasing  proceeds  as 
follows: 

1)  If  Pvf  >  10  P!,,  velocities  are  from  first  trip  clutter  residue 

and  first  trip  echo.  (Calculation  of  reflectivity 

should  be  from  velocity  channel.) 

2)  If  2PJ,  <  pyf  <  lOP^,  recovery  of  velocities  is  not  possible.  (Reflec¬ 

tivities  calculated  from  reflectivity  channel.) 

3)  If  P^  <  2P'^  velocities  are  from  the  second  trip  signal.  (Reflec¬ 

tivities  calculated  from  velocity  channel.) 

Note  that  in  this  interlaced  mode,  assignment  of  ranges  to  the  first  trip  is 
very  satisfying  since  errors  of  the  order  of  1  dB  in  P^  and  P£  do  not  signifi¬ 
cantly  affect  the  comparison  #1.  Not  so  in  the  comparison  #3.  Ideally,  one  would 
want  to  have  P^  plus  the  clutter  residue  power  Pr  smaller  by  10  dB  from  P^.  But 
because  the  filter  gain  in  the  passband  varies  by  ±1  dB  and  because  of  previously 
mentioned  uncertainties,  comparison  #3  is  chosen.  Therefore,  when  condition  #3  is 
satisfied,  there  will  be  times  when  erroneous  velocities  are  obtained  because 
Pr+P-j  may  not  be  smaller  than  0.1  P2.  On  the  other  hand,  if  the  condition  is 
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and  gain  is  near  unity. 


REFLECTIVITY  CHANNEL 


Figure  7.  Block  diagram  of  the  reflectivity  and  velocity  channels 
with  ground  clutter  filters. 


changed  to  Pf  <  1.1  P^.  more  data  will  be  lost,  and  still  some  may  be  erroneous 
due  to  the  statistical  and  systematic  uncertainties.  A  phase  diversity  radar 
would  help  in  this  instance  because  biases  (due  to  larger  values  of  Pr+P1  allowed 
by  #3)  would  be  eliminated  at  the  expense  of  increased  standard  errors  of  estimates. 

Our  discussion  carries  to  other  choices  of  the  two  PRT's.  For  instance,  if 
reflectivity  is  estimated  with  a  PRT  of  3TS>  one  needs  to  make  a  comparison  in  the 
reflectivity  channel  of  powers  in  the  second  P£  and  third  P^  ambiguous  range  ring 
for  velocity  estimation.  So  wherever  P£  £  10  P^,  assume  significant  power  is  in 
second  trip,  and  compare  P,!,+Pg  with  P  ^  as  in  #1  to  #3.  Whenever  P^  >  10  P£, 
assume  significant  power  is  in  third  trip,  take  P^+P^  and  proceed  with  #1  to  #3. 

If  P 2  and  P^  are  within  10  dB  of  each  other,  data  in  second  and  third  trip  are 
contaminated,  take  P^+P^  and  compare  as  in  #1  to  #3. 

Although  recovery  of  the  second  trip  storm  velocities  if  overlaid  by  ground 
clutter  and  first  trip  storm  echo  is  not  satisfactory,  the  overall  performance  of 
the  radar  will  be  good.  First,  we  must  bear  in  mind  that  occurrence  of  triple 
overlay  is  rare.  Second,  when  only  ground  clutter  overlays  the  storm  in  the 
second  trip,  a  good  canceler  (50  dB)  will  reduce  clutter  residue  to  a  negligible 
value  in  most  of  the  area  so  that  comparison  #3  will  produce  valid  data. 

An  alternate  method  to  assign  velocities  is  by  comparison  of  powers  in  the 
reflectivity  channel  after  the  filter.  This  is  equivalent  to  the  present  scheme 
on  NSSL  radars  (except  that  they  don't  yet  have  a  clutter  filter  in  the  reflec¬ 
tivity  channel).  If  such  a  comparison  were  made,  the  full  benefit  of  the  filter 
in  the  velocity  channel  would  not  be  real i zed; a  clutter  filter  in  the  reflectivity 
channel  (when  interlaced  samples  are  used)  would  cancel  less  than  the  filter  in  the 
velocity  channel  because  the  number  of  samples  for  reflectivity  measurements  is 
several  times  smaller  than  for  velocity  measurements. 

4.  The  Third  Order  Recursive  Filter 

For  accurate  mean  velocity  estimation  (rms  error  51  m*s"^)  via  the  pulse  pair 
processor,  one  needs  a  10  dB  ratio  of  signal  to  clutter  residue  and  15  dB  ratio  if 
spectrum  width  is  estimated  (Zrnic1  and  Hamidi,  1981).  Thus,  the  canceler  must 
reduce  the  residue  to  at  least  10  dB  below  signal  power.  Therfore,  a  50  dB 
canceler  would  allow  moment  estimation  up  to  about  40  dB  of  clutter  C  to  signal  S 
power  ratio. 

We  consider  a  50  dB  rejection  in  a  stopband,  a  1  dB  ripple  in  the  passband, 
and  a  transmission  zero  at  DC,  reasonable  goals  to  achieve  since  a  recursive  filter 
with  such  characteristics  has  already  been  designed  and  fabricated  by  Raytheon 


engineers  (Groginsky  and  Glover,  1980).  Incidentally,  no  more  than  50  dB  of  can¬ 
cellation  is  possible  if  an  8-bit  A/D  converter  is  used  because  when  the  clutter 
spans  the  full  8-bit  range,  its  quantization  noise  is  50  dB  below  its  average  power. 

In  Table  3,  we  list  a  third  order  filter  that  we  have  analyzed,  and  on  Figure 
8  is  the  appropriate  flow  chart  that  defines  the  k  coefficients.  For  illustration, 
the  transfer  characteristic  of  the  filter  is  shown  in  Figure  9  where  the  stopband 
cutoff  velocity  v$  and  passband  cutoff  velocity  vp  are  also  defined. 

The  main  disadvantage  of  the  recursive  elliptic  filter  is  its  transient 
response.  Although  the  steady  state  frequency  response  of  HR  (Infinite  Impulse 
Response)  filters  may  compare  well  with  FIR  (Finite  Impulse  Response)  filters, 
infinite  duration  of  their  impulse  response  has  limited  use  in  radar  application. 

At  low  elevation  angles  (0.5°)  the  clutter  power  is  several  orders  of  magnitude 
stronger  than  the  signal  power.  Thus,  when  interlaced  samples  (batch)  for  velocity 
and  reflectivity  estimation  are  used  (Doviak  et  al . ,  1979),  the  clutter  appears  as 
a  step  to  the  filter  which  produces  severe  ringing  at  the  output  and  in  effect 
overwhelms  the  signal.  A  similar  condition  occurs  if  the  automatic  gain  control 
(AGC)  changes  the  gain  just  before  processing  a  block  of  pulses.  To  reduce  this 
effect,  either  enough  time  must  be  allowed  for  the  transients  to  settle  before 
pulse  pair  processing  is  attempted  or  some  other  means  of  their  suppression  must  be 
employed.  Because  settling  time  is  tens  of  pulses,  we  consider  the  alternative 
which  is  initialization  (Fletcher  and  Burlage,  1972).  The  purpose  of  initializa¬ 
tion  is  to  load  the  memory  elements  of  the  recursive  filter  with  their  anticipated 
steady  state  values.  These  values  can  be  estimated  from  the  first  incoming  pulse. 
When  the  ground  clutter  is  a  perfect  D.C.  line  and  no  weather  signal  is  present, 
the  initialization  will  immediately  bring  the  filter  into  steady  state.  However, 
the  clutter  has  finite  width  and  superposed  on  it  is  the  weather  signal  which 
together  degrade  the  initialization. 

For  the  third  order  filter  on  Figure  8,  the  output  of  the  first  memory  element, 
Z-1,  would  be  set  to  A^l-K^)  where  A^  is  the  digital  value  of  the  first  sample. 

The  outputs  of  the  other  two  elements  must  be  set  to  0.  With  this,  the  output 
after  the  first  time  sample,  A-j,  is  0  and  after  the  second  one,  A2»  it  is  A2~A^ . 
Thus,  the  filter  immediately  acts  as  a  delay  line  canceler  so  that  a  large  D.C. 
value  does  not  get  through. 

The  expected  clutter  attenuation  for  an  initialized  filter  is  25  dB  more  (in 
transient)  than  for  the  noninitial ized  filter  (Figure  10).  Attenuations  were 
computed  from  the  expected  power  at  the  filter's  output  for  a  given  recursion. 
Increase  in  spectrum  width  reduces  the  attenuation  because  more  skirts  of  the 
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Table  3. 


Coefficients  for  the  high  pass  elliptic  filter  (PRT=781  ys,  X  =  10  cm). 


Hz  (m/s)  Hz  (m/s) 


60  2.9  17.0  0.8  1.995322175  1.802510637  0.871228116  0.577292360 


clutter  spectrum  spill  into  the  passoand.  It  is  noteworthy  that  the  transient  for 
the  initialized  filter  lasts  about  50  recursions  but  for  the  noninitialized  it  is 
over  100.  Moreover,  the  initialized  filter  can  be  quite  effective  even  in  tran¬ 
sient  because  the  attenuation  is  considerable.  At  the  first  recursion  the  non¬ 
initialized  filter  does  not  attenuate  clutter  due  to  a  direct  path  between  input 
and  output  (Figure  8).  This  is  not  so  for  the  initialized  filter  where  the  output 
starts  at  a  second  recursion. 

In  the  interim  report  (Zrnic  and  Hamidi,  1981),  the  following  equation 
relating  the  clutter  spectrum  width  oc,  the  passband  velocity  vp  and  the  frequency 
domain  rejection  L  (dB)  was  derived: 

vp  ~  °c  5  log  e  *  ^ 

For  a  50  dB  rejection,  this  becomes 

vp  >  4.8  ac  (3) 

Using  calculations  that  lead  to  Figure  10,  we  have  obtained  anothet  rela¬ 
tionship  between  a  and  v  for  a  fixed  integrated  attenuation.  The  attenuation  in 
c  p 

steady  state  decreases  with  the  increase  in  clutter  spectrum  width  (Figure  11), 
and  we  note  50  dB  or  more  attenuation  is  for  widths 

°c  <  vs  (4) 

For  this  type  of  filter  vp  =  3.6  v$  and  in  order  to  keep  the  total  residual 
power  below  50  dB  we  find 

*p  >  8-2  °c  (5) 
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ATTENUATION  (dB) 


Figure  8.  Block  diagram  of  the  third  order  high  pass  elliptic  filter 


VELOCITY  (m  r') 


Figure  9.  Frequency  (Doppler  velocity)  response  characteristic  of 
the  recursive  filter.  The  Byquist  interval  is  64  m- s~*. 
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Figure  10.  Expected  attenuation  after  a  clutter  step  is  applied  to  the 
filter.  Two  clutter  spectrum  widths t  aa  =  0.2  m-s~ and  aa  =  0.5  nrs~l 
are  considered  for  both  the  initialized  and  noninitialized  filter. 


The  recursive  filter  operates  on  each  component  of  the  complex  digital  video 
samples  before  the  pulse  pair  processor  (Figure  12).  In  addition  to  the  mean 
velocity  v  and  spectrum  width  a  estimates,  the  processor  provides  an  estimate  of 

a  * 

echo  power  P.  Power  is  computed  from  a  uniform  average  of  the  square  of  time 
sample  magnitudes.  Mean  velocity  is  from  the  argument  of  the  autocovariance  at 
lag  one  and  the  spectrum  width  from  the  logarithm  of  the  correlation  coefficient 
at  lag  one.  For  short,  we  call  all  three  pulse  pair  processing. 

To  illustrate  the  performance  of  the  initialized  tandem  filter-pulse  pair 
processor,  simulations  on  synthetic  time  series  data  and  on  some  real  data  were 
conducted.  Roundoff  effects  on  coefficients  and  finite  word  lengths  of  data  were 
not  considered.  Groginsky  and  Glover  (1980)  show  that  coefficients  with  12-  to 


SPECTRUM  WIDTH  <rc  (m  »-') 


NORMALIZED  CLUTTER  SPECTRUM  WIDTH  <rc\ 

Figure  11.  Expected  clutter  attenuation  for  a  filter  in  steady  state. 


Figure  12.  Block  diagram  of  data  flow.  I  and  Q  are  the  digitized  video 
signals  that  must  he  filtered  prior  to  pulse  pair  processing. 


16-bit  word  lengths  (depending  on  radar  PRT)  do  not  degrade  the  filter  performance 
(see  section  5). 

Digital  computer  simulation  of  the  block  diagram  on  Figure  12  tested  the 
quality  of  the  output  after  filtering.  The  input  time  series  data  were  synthe¬ 
sized  such  that  they  contained  a  clutter  and  weather  signal  of  Gaussian  spectral 
shapes  and  had  Gaussian  statistics  (see  Zrnic',  1975,  for  details).  Unambiguous 

velocity,  v  .  was  set  to  32  m-s"1. 
a 

In  the  results  that  follow,  eight  complex  time  samples  were  passed  through 

A  A  A 

the  filter  and  estimates  P,  v,  and  were  computed.  We  consider  eight  contiguous 
pulses  to  be  close  to  a  lower  limit  as  far  as  standard  errors  of  estimates  is 
concerned.  These  errors  are  less  than  about  3  m-s”1  in  both  spectrum  width  and 
velocity  if  widths  are  less  than  8  m-s"1  (Zrnic1,  1979b).  The  actual  number  of 
contiguous  pulse  pairs  out  of  eight  pulses  becomes  six  since  the  first  pulse  is 
lost  in  the  process  of  initialization. 

Even  if  3  m-s  1  errors  are  not  tolerable,  groups  of  eight  pulses  might  be 
used  in  an  interlaced  pulse  scheme  to  estimate  velocity  and  spectrum  width.  The 
reflectivity  is  obtained  from  a  separate  pulse  with  a  clearing  period  to  achieve 
large  unambiguous  range.  Comparison  of  powers  in  the  reflectivity  channel  enables 
one  to  assign  correct  ranges  to  the  velocity  estimates  that  have  a  small  unambig¬ 
uous  range  (Doviak  et  al . ,  1979).  This  technique  requires  that  coherent  and  ’■ 
coherent  estimates  are  from  echoes  that  originate  in  overlapping  resolution 
volumes;  a  small  number  of  contiguous  samples  for  velocity  and  -pectnM*;  width 
estimates  insures  enough  overlap.  Hence,  our  choice  of  ei gr*  pulses  is  justified 
on  one  more  ground. 

A  A 

Mean  values  of  the  estimates  P,  v,  and  oy  can  serve  to  quantify  the  compari¬ 
son  between  a  noninitial ized  and  initilized  filter  because  these  estimates  are  the 
final  product  of  interest.  Blocks  of  8  time  samples  were  applied  repeatedly  100 
times  to  the  recursive  filter  and  that  many  estimates  of  power,  velocity  and 
spectrum  width  were  obtained.  The  estimates  were  then  averaged  to  generate  the 

mean  values.  Throughout  this  simulation  three  clutter  widths,  a  =  0.2,  0.4,  and 
-1  _  I  c 

0.6  m-s  were  examined  because  the  first  one,  0.2  m-s ",  is  much  narrower  than 

the  filter  notch;  the  second,  0.4  m-s"1,  is  only  slightly  larger  (by  0.05  m-s"1) 

than  the  value  allowed  by  inequality  (3);  and  0.6  m-s”1  satisfies  the  inequality 

in  (5).  Also,  the  percentage  of  occurrences  of  clutter  widths  below  0.6  m-s"1 

(for  the  Norman  radar)  is  about  constant,  but  drops  rapidly  after  0.6  m-s"1  if 

the  antenna  is  rotating  at  10°s  (see  Figures  1.9,  1.10,  and  1.11  in  Zrnic'  and 

Hamidi,  1981). 
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Cl  utter- to- signal  ratios  of  10,  20,  30,  and  40  dB  for  the  initialized  filter, 
and  ratios  of  -10,  0,  10,  and  20  dB  for  the  noninitialized  one  were  examined.  As 
expected,  power  estimates  from  the  initialized  filter  are  less  than  one  dB  from 
their  true  values  as  long  as  C/S  -  30  dB  and  the  mean  velocity  is  larger  than 
8  m>s”^  (Figure  13).  At  a  narrow  signal  spectrum  width  (e.g.,  1  m-s~^),  near  zero 
velocity  the  filter  removes  quite  a  bit  of  power,  but  at  a  width  of  8  m*s"\  at 
most  3  dB  is  lost.  Clutter  power  must  be  reduced  by  about  20  dB  in  order  to 
achieve  comparable  performance  with  a  noninitialized  filter. 

Mean  velocity  estimates  (Figure  14)  also  illustrate  the  improved  performance 
of  the  initialized  filter.  The  first  sample  was  dropped  and  six  pairs  (out  of 
seven  possible)  were  used  in  order  to  have  the  same  number  for  two  filter  condi¬ 
tions  (i.e.,  initialized  and  noninitialized).  Note  that  the  initialized  filter  is 
effective  up  to  20  dB  of  cl  utter- to-signal  ratio,  which  is  20  dB  better  than  the 
noninitialized  case.  A  slight  positive  bias  at  small  velocities  and  weak  clutter 
is  due  to  loss  of  the  portion  of  weather  spectrum  near  the  origin.  Otherwise, 
strong  clutter,  30  or  more  dB  above  signal,  creates  unacceptable  bias.  Thus,  as 
far  as  mean  velocity  bias  is  concerned,  this  filter  with  a  50  dB  notch  when 
properly  initialized  acts  like  a  30  dB  filter  (which  includes  the  10  dB  margin 
discussed  in  the  beginning  of  this  section).  Without  initialization,  the  effective 
attenuation  is  only  10  dB,  which  can  be  deduced  from  the  fact  that  at  C/S  =  0  dB 
a  good  mean  velocity  estimate  (bias  about  1  m*s_1)  is  obtained.  To  achieve  similar 
results  without  the  filter,  clutter  power  must  be  10  dB  below  signal. 

Spectrum  width  estimates  had  unacceptable  bias  that  depends  on  clutter-to- 
signal  ratio,  mean  velocity,  and  spectrum  width.  For  instance,  with  mean 
velocity  of  24  m-s-1  (i.e.,  3/4  Nyquist  velocity)  and  C/S  *  10  dB,  bias  was  larger 
than  3  m*s”^.  Such  large  bias  is  due  to  the  transients  (damped  sinusoid)  of  the 
second  order  section.  These  transients  create  symmetrical  sidebands  around  zero 
velocity. 

It  is  possible  to  reduce  the  disastrous  effect  of  transient  on  spectrum  width 
by  having  a  longer  train  of  pulses  and  ignoring  several  of  those  at  the  beginning. 
Figure  10  may  guide  the  determination  of  the  number  of  samples  that  should  be 
allowed  for  transients  to  settle,  In  order  to  obtain  good  clutter  canceling.  For 
Instance,  if  the  first  16  filtered  samples  are  dropped,  the  clutter  with  a  spec¬ 
trum  width  of  0.2  m*s”^  is  attenuated  by  over  36  dB.  The  examples  on  Figure  15  was 
obtained  when  32  pulses  were  passed  through  the  initialized  filter  and  the  last 
sixteen  are  used  for  width  calculation.  Results  are  satisfactory  (biases  less 
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Figure  IS.  Power  estimate  after  the  noninitiaiized  recursive  filter  (upper  graphs) 
and  initialized  filter  (lower  graphs).  On  the  ordinates  are  logarithms  of  P/P 
where  P  is  the  true  mean  power.  Clutter  width  CLWD  is  0. 4  m •  signal  widths 

SIGWD  are  1  m*  s~^  (on  the  left)j  and  8  m’8~ ^  (on  the  right) >  and  clutter- to- signal 
power  ratio  GC/SIG  is  indicated.  Eight  samples  were  passed  through  the  filter , 
hut  only  the  last  seven  were  used  for  power  estimation. 
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Figure  14.  Estimated  mean  velocity  versus  true  mean  velocity  for  simulated  times 
series  data  processed  by  the  pulse  pair  processor.  Other  parameters  are  as  on 
Figure  13. 
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Figure  15.  Estimated  mean  signal  spectrum  width  versus  true  width  for  simulated  time  series  data.  Thirty-two 
samples  were  filtered  through  the  initialized  filter  hut  only  the  last  sixteen  were  applied  to  the  pulse  pair 
processor.  Clutter  widths  are  0.2  m’S~l  (upper  graphs)  and  0.4  m-  s~^  (lower  graphs).  The  mean  signal 
velocities  from  left  to  right  are  2 ,  8}  and  24  m*  s--2. 


than  1  m»s~  )  for  cl  utter- to-signal  ratios  of  10  and  20  dB,  but  fall  apart  at  30  dB 
and  mean  velocities  larger  than  12  m-s"1. 

5.  Design  Consideration 

The  representation  of  filter  coefficients  and  the  architecture  of  the  third 
order  recursive  filter  are  briefly  discussed  in  this  section. 

Filter  coefficients  and  data  in  our  implementation  will  be  converted  to 
fixed-point  fractional  numbers.  Figure  16  shows  the  response  of  the  filter  with 
coefficients  represented  as  fixed-point  numbers  with  10-bit  fractional  parts, 
along  with  the  response  of  the  actual  filter.  As  can  be  seen,  the  filter  per¬ 
formance  is  quite  poor  since  only  30  dB  of  cancellation  is  achieved.  Th*  perfor¬ 
mance  at  wider  stopbands  is  somewhat  better.  Figure  17  shows  the  response  of  a 
filter  with  12  bits  in  fractional  part,  which  meets  the  design  objectives. 

Figure  18  is  the  block  diagram  of  the  proposed  hardware  implementation  of  the 
three-pole  high  pass  elliptic  filter  used  for  ground  clutter  cancellation.  The 
hardware  is  based  on  a  high-speed  parallel  arithmetic  and  logic  unit  assisted  by 
four  multipliers  capable  of  simultaneous  operation.  The  hardware  is  capable  of 
performing  the  entire  filtering  arithmetic  within  1  microsecond,  resulting  in  real¬ 
time  operation.  System  control  is  implemented  in  microcode  to  maintain  flexibility 
and  facilitate  future  modifications  if  required.  The  control  logic  operates  at  a 
clock  rate  of  10  MHz.  Filter  coefficients  and  constants  are  implemented  in 
switches,  and  the  filtering  range  can  also  be  selected  by  switches.  Filter 
initialization  is  provided  in  hardware  and  is  automatically  performed  when  the 
radar  is  operating  in  the  expanded  (interlaced  PRF)  mode. 

The  hardware  is  implemented  with  a  16-bit  parallel  internal  data  bus.  Data 
and  coefficients  are  represented  by  fixed-point  binary  numbers  with  12-bit  frac¬ 
tional  part. 

The  design  and  fabrication  of  the  filter  are  being  supported  by  Joint  Systems 
Program  Office,  and  we  anticipate  that  the  device  will  be  operational  in  the  first 
half  of  1983. 

6.  Conclusions 

Problems  associated  with  the  first  trip  clutter  overlaying  second  trip  area 
for  velocity  and  spectrum  width  estimation  have  been  pointed  out.  The  obscured 
area  in  the  second  trip  is  larger  (by  about  10  times  for  the  NSSL  radar)  than  the 
first  trip  clutter  area.  Moreover,  the  weather  signals  in  the  second  trip  area 
suffer  a  range  square  disadvantage  which  for  the  NSSL  radar  is  about  10  to  25  dB. 
The  very  important  aspect  concerning  range  de-aliasing  of  velocity  estimates  has 
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Figure  16.  Frequency  (velocity)  response  of  a  narrow  notch  filter.  Solid  line  is 
for  ideal  coefficients  and  dashed  line  is  with  a  10-bit  representation. 
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Figure  17.  Frequency  response  of  the  filter  on  Figure  16  except  a  12-bit  repre 
sentation  is  used. 


been  examined.  Much  work  remains  to  be  done  in  this  area  to  realize  the  full 
benefit  of  the  ground  clutter  canceler  in  the  velocity  channel.  The  optimum,  but 
expensive,  solution  calls  for  a  dual  frequency  radar  with  two  coherent  receivers 
and  ground  clutter  filter  with  almost  equal  characteristics.  Then  the  powers  in 
the  reflectivity  channel  can  be  used  to  assign  correct  ranges  for  the  velocity 
estimates.  But,  with  interlaced  sampling,  clutter  filters  in  the  reflectivity  and 
velocity  channel  would  have  quite  different  characteristics;  the  notch  depth  would 
be  probably  20  dB  deeper  in  the  velocity  channel  because  the  dwell  time  for 
velocity  estimation  is  longer.  Consequently,  when  powers  in  the  reflectivity 
channel  are  compared  for  assigning  correct  ranges  to  the  velocities,  the  full 
potential  of  the  clutter  filter  in  the  velocity  channel  will  not  materialize,  i.e., 
some  second  trip  velocities  would  be  falsely  declared  invalid.  A  more  complicated 
processing  requires  comparisons  of  powers  in  the  reflectivity  channel  with  those 
in  the  velocity  channel. 

The  most  economical  canceling  scheme  is  with  a  third  order  recursive  filter. 
This  filter  achieves  a  50  dB  rejection  in  the  stopband  with  total  annihilation  of 
DC.  One  dB  ripple  is  in  the  passband,  and  the  ratio  of  passband  cutoff  vp  to 
stopband  cutoff  v$  is  about  3.6.  The  filter  operates  best  in  steady  state,  but  it 
can  also  be  made  to  operate  in  transient  by  properly  initializing  its  memory 
elements.  Performance  on  8  simulated  time  samples  shows  that  about  10  dB  of 
cl  utter- to- signal  margin  is  lost  for  power  and  20  dB  for  mean  velocity  estimation. 
But  with  such  a  small  number  of  samples  the  width  estimate  is  useless.  Thus, 
interlaced  sampling  can  be  used  for  velocity  estimation  at  a  20  dB  penalty  in 
clutter  canceling,  but  it  is  not  suitable  for  spectrum  width  estimation. 

In  order  for  spectrum  width  biases  to  be  less  than  1  m*s_1,  a  longer  train  of 
pulses  must  be  employed  on  an  initialized  filter  and  several  of  the  leading  pulses 
must  be  used  exclusively  to  further  reduce  the  transients.  Then  from  the  rest  of 
the  pulses  the  widths  will  be  credible.  Theoretical  curves  of  the  expected  atten¬ 
uation  for  the  initialized  filter  in  transient  are  presented.  From  these  one  can 
determine  the  number  of  recursions  needed  to  reduce  clutter  power  below  a  desired 
level  and  thus  the  number  of  samples  that  should  not  be  applied  to  the  pulse  pair 
processor.  An  example  shows  that  if  the  first  16  samples  out  of  32  are  solely 
used  to  further  reduce  the  transients,  spectrum  width  can  be  estimated  with  little 
bias  as  long  as  clutter- to-signal  ratio  is  less  than  20  dB. 

The  part  of  weather  signal  that  is  within  the  stopband  cutoff  of  the  filter 
is  totally  lost,  but  that  is  a  small  penalty  because  zero  radial  velocities  occupy 
a  small  portion  of  a  storm  and  should  not  affect  significantly  total  areal 


rainfall  rate  estimates.  Furthermore,  in  these  instances  the  mean  velocities  will 
be  noisy  but  not  biased.  However,  large  spectrum  widths  indicative  of  turbulence 
and  other  hazards  would  still  be  measurable.  Therefore,  canceling  should  be  done 
selectively  at  those  ranges  and  elevation  angles  where  ground  clutter  is 
significant. 

Examination  of  hardware  for  implementing  the  filter  resulted  in  a  preliminary 
design  based  on  a  high  speed  parallel  arithmetic  and  logic  unit  assisted  by  four 
multipliers  capable  of  simultaneous  operation.  Two  identical  units  are  needed, 
one  for  the  inphase,  the  other  for  the  quadrature  component  of  the  video  signal. 
The  proposed  architecture  has  a  throughput  rate  faster  than  1  ps  and  is  thus 
capable  of  real  time  operations. 

We  stress  that  clutter  canceling  is  easiest  with  an  uninterrupted  pulse  train 
because  the  recursive  filter  is  then  less  subject  to  transients. 
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